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We study the quantum dynamics of a Michelson interferometer with Fabry-Perot cavity arms and one movable 
end mirror, and driven by a single photon — an optomechanical device previously studied by Marshall et al. 
as a device that searches for gravity decoherence. We obtain an exact analytical solution for the system's 
quantum mechanical equations of motion, including details about the exchange of the single photon between 
the cavity mode and the external continuum. The resulting time evolution of the interferometer's fringe visibility 
displays interesting new features when the incoming photon's frequency uncertainty is narrower or comparable 
to the cavity's line width — only in the limiting case of much broader-band photon does the result return to 
that of Marshall et al, but in this case the photon is not very likely to enter the cavity and interact with the 
mirror, making the experiment less efficient and more susceptible to imperfections. In addition, we show that 
in the strong-coupling regime, by engineering the incoming photon's wave function, it is possible to prepare the 
movable mirror into an arbitrary quantum state of a multi-dimensional Hilbert space. 

PACS numbers: 



I. INTRODUCTION 

Recently, significant progress has been made in observ- 
ing quantum effects in macroscopic mechanical systems [1]. 
As presented in the work of O'Connell et al. [2], a 6-GHz 
nano-mechanical oscillator was cooled down near its quantum 
ground state with dilution refrigeration, and later prepared 
into a Fock state by coupling the oscillator to a supercon- 
ducting qubit. States with thermal occupation numbers below 
unity have also been achieved with cavity-assisted radiation- 
pressure cooling, by Teufel et al. [3] and Safavi-Naeini et 
al. [4]. Further more, as shown by Gupta et al. [5] and Thomp- 
son et al. [6], it is possible to couple a single photon strongly 
with a mechanical degree of freedom, such that the momen- 
tum imparted by a single photon to a mechanical degree of 
freedom can be comparable to its initial momentum uncer- 
tainty. 

In this paper, we study the open quantum dynamics of a 
nonlinear optomechanical device, namely a Michelson inter- 
ferometer with Fabry-Perot cavities, one of them with a mov- 
able end mirror (acting as the mechanical oscillator). This 
device, driven by a single photon, was proposed by Marshall 
et al. [7, 8] as an experiment to search for Penrose's conjec- 
ture of gravity decoherence [9]. Such single-photon driven 
devices have also been more recently studied by Rabl [10] 
and Nunnenkamp et al. [11]. By taking advantage of the con- 
served quantity — the total number of photons in the system, 
one can obtain exact solutions to this system's quantum dy- 
namics. Unlike Rabl and Nunnenkamp et al., who studied 
systematically the statistics of the out-going photons and the 
steady state reached by the mechanical oscillator, we focus in- 
stead on the fringe visibility of a single-photon interferometer, 
and the conditional quantum state of the mechanical oscillator 
upon the detection of an out-going photon. 

The single-photon Michelson interferometer is shown 
schematically in Fig. 1, in which the port on the left is the 
input port, towards which the single photon is injected; the 
photon, after interacting with the Michelson interferometer, 
may exit either from the input port, or from the other open 
port. Each of the two arms consists of a high-finesse opti- 



cal cavity; the setup of these two cavities are identical, ex- 
cept one of them has a movable end mirror, which acts as 
the mechanical oscillator that interacts with light in the cavity. 
The 50/50 beam splitter splits the quantum state of the entire 
mirror-light system into two components, one of them corre- 
sponding to the photon entering the fixed cavity (and leaveing 
the oscillator at its initial state), the other corresponding to the 
photon entering the movable cavity (thereby modifying the 
oscillator's state through radiation pressure). We will set the 
displacement zero-point of the interferometer to have equal 
arm lengths, with each arm at a distance equal to the beam- 
splitter. At such a zero point, the photon injected from the 
input port will return to the input port with unit probability. 
Therefore we also call the input port the "bright port" and the 
other open port the "dark port". We can artificially tune the 
interferometer away from its zero point, e.g., by adjusting the 
fixed microscopic distances between the front mirrors and the 
beamsplitter. This changes the relative phase <p between the 
two superimposed components in wave function of the entire 
system; the resulting variations in the probability density of 
having the photon exiting the bright port at time t, quantified 
by the fringe visibility, is a measure of the degree of coherence 
between these two components at this moment in time. 

In the case of low environmental temperature and in the ab- 
sence of unexpected mechanisms of decoherence, Marshall et 
al. showed that the visibility will revive completely for every 
half of the mechanical oscillation period. In obtaining such 
a result, they assumed the photon was initially already in ei- 
ther of the two cavity arms, and considered a closed evolu- 
tion of the cavity mode and the mechanical oscillator. This 
assumption has also been widely used in analysis of such a 
nonlinear optomechanical device, e.g., by Bose et al. [12] and 
subsequent analysis of the Marshall experiment [7] by Bassi 
et al [13]. 

In a realistic experimental setup, it is necessary to take full 
account of the open quantum dynamics of this system, which 
that involves the oscillator (the mirror), the cavity mode and 
the external continuous field, including how the single photon 
is coupled into the cavity in the first place. The open quantum 
dynamics depends on the wave function of the photon, whose 
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Figure 2: (color online) A schematics showing a single-photon in- 
terferometer with Fabry-Perot cavity and a movable mirror. The dis- 
placement of the mirror-endowed mechanical oscillator y is paramet- 
rically coupled to the cavity mode a, which has an eigenfrequency 
C0q with y = 0. The cavity mode in turn couples to the ingoing con- 
tinuous filed c(x) and outgoing continuous field d(x). 



Figure 1 : (color online) A schematics showing the single-photon in- 
terferometer. The external single photon excites the cavity mode 
which in turn interacts with the movable end mirror via radiation 
pressure. This is adapted from Fig. 1 of Ref. [7] with small modifi- 
cations. 



coupling between the mirror and the cavity mode is paramet- 
ric: the position y of the mirror modifies the eigenfreqeuncy 
of the cavity mode. 



Fourier transform is related to the frequency content of the 
photon. For example, if the photon has a short-pulse wave 
function with time-domain duration much less than the cavity 
storage time, which corresponds to a frequency uncertainty 
much larger than the cavity line width, then the photon will 
only enter the cavity with a small probability. By contrast, a 
narrowband photon (with frequency uncertainty below cavity 
line width) must have a wave packet duration much longer 
than cavity storage time, and therefore we must address the 
issue that the photon can be simultaneously inside and outside 
the cavity. The latter scenario, although more complicated, 
might be experimentally more favorable, as in this scenario 
the photon has a high probability to enter the cavity and to 
interact with the mirror much more strongly. 

The outline of this article goes as follows: in Sec. 11, we will 
write down the Hamiltonian of our nonlinear optomechanical 
device and study the open quantum dynamics by solving the 
Shrodinger equation exactly; in Sec. Ill, we will give a de- 
tailed analysis of the single-photon interferometer, and will 
calculate the interferometer's fringe visibility; in Sec. IV, we 
will show that the mechanical oscillator can be prepared to an 
arbitrary quantum state in a multi-dimensional Hilbert space, 
if we inject the single photon with a properly-designed pro- 
file into the interferometer; in Sec. VI, we will summarize our 
main results. 



II. A SINGLE CAVITY WITH ONE MOVABLE MIRROR 

Before studying the entire single-photon interferometer, we 
first consider a single cavity, as shown schematically in Fig. 2. 
The cavity has one fixed mirror located at x = 0, and one mov- 
able mirror which acts as a mechanical oscillator. Here, as- 
suming the injected photon to have a frequency content much 
less than the free spectral range of the cavity (which has a rel- 
atively high finesse), we will only consider one optical mode 
of the cavity (which we shall refer to as the cavity mode). By 
assuming a high finesse for the cavity, this mode couples to 
the external vacuum via single-photon exchange. At linear 
order in the mirror's motion and assuming low velocity, the 



A. The Hamiltonian 

Here we will write down the Hamiltonian of the system. 
For simplicity, we will use natural units with h = 1 and c = 1 
throughout this paper. The Hamiltonian of the external con- 
tinuous optical field, in the position space representation, is 
given by 

[(d x c T x )c x ~c T x d x c x ]dx 



H, 



{(d x dl)d x - d\d x d x ] dx 



(1) 



where c x and d x are the annihilation operators for ingoing and 
outgoing field at location x, respectively. Note that for the 
actual setup shown in Fig. 2, the ingoing and outgoing field 
are on the same side of the front mirror, namely both at x < 0. 
Since the field operators at different locations commute with 
each other — [c x cj,] — 8(x — x'), we can fold the outgoing field 
from [— °o,0] into [0,+oo], therefore just use c to denote both 
the ingoing and outgoing fields, with c x (x < 0) for the ingoing 
field and c x (x > 0) for the outgoing field, namely 



Hn 



(d x clc x -cld x c x )dx. 



(2) 



The free Hamiltonian of the single cavity mode is given by 



H c = (£>Qa* a. 



(3) 



with a the annihilation operator and [a, a'] = 1. 

The free Hamiltonian for the mechanical oscillator reads 



H„, 



Py , 1 „2 -2 
2m + 2 m( °^ 



(4) 



where y and p y are the position and momentum operators, re- 
spectively. 

The total interaction Hamiltonian Hi between the external 
continuum and the cavity mode in the rotating-wave approxi- 
mation, and between the cavity mode and the mechanical os- 
cillator, is given by 



Hi = iy/Y(codf — acj) + led' ay. 



(5) 
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Here y = £■ is the cavity bandwidth with L being the cav- 
ity length; k = (Oq/L is the optomechanical coupling constant. 
The interaction between the cavity mode and the external con- 
tinuum takes place at the front mirror with x = and the 
Hamiltonian describes the exchange of photon between them. 
The total Hamiltonian is a sum of the free and the interaction 
parts, namely, 

H = H +H c +H m +H,. (6) 

Note once more that by including only a single cavity mode 
resonant at frequency (Oq/(2tz), we must make sure the fre- 
quency content of the injected light is focused well within a 
free spectral range, c/(2L). 

B. Structure of the Hilbert Space 

Even though the Hamiltonian contains a cubic term a* ay, 
which implies a nonlinear dynamics, we have a conserved dy- 
namical quantity — the total photon number: 

d'd + J c x c x dx, (7) 

which makes the system's evolution still analytically solvable, 
as also recognized by Rabl [10] and Nunnenkamp et al. [11]. 
Since the initial state of our system consists of one single pho- 
ton, there can only be one photon throughout the entire evo- 
lution. Mathematically, this means we only need to consider 
a one-photon subspace of the entire Hilbert space, which in 
turn consists of three disjoint subspaces, which corresponds 
to: J%i-, which corresponds to an incoming photon towards 
the cavity; J%, which corresponds to a photon inside the cav- 
ity, and ,3#\+, which corresponds to a photon leaving the cav- 
ity. All quantum states in this space can be written as: 

W) = 'f(x,t)e- i<aa <'- x )\x) 7 ®\<l> l (x,t)) m dx 

+ O(0« _ft * , a t |0) r (»|fc(0)m. (8) 

Here 

\x) 7 =cl\0) 7 (9) 

is the "position eigenstate" of the single photon outside of 
the cavity, and |0)y is the optical vacuum; the subscripts 
7 and m indicate Hilbert spaces of light and movable mir- 
ror, respectively; f(x,t) is a complex function of position 
(— °° < x < +°°) and time, a(t) is a complex function of time t 
alone; |0i (x,t)) m and |02(O)m are two families o/state vectors 
that belong to the Hilbert space of the mechanical oscillator. 
At any given time, the x < part of the integral term on the 
right-hand side corresponds to J%\ _ , the x > part of the in- 
tegral term corresponds to while the non-integral term 
corresponds to Ji%. In general, all three terms will be present, 
which means the entire system's quantum state is a superpo- 
sition of having the photon simultaneously present in all three 
possible locations. Note that the factors e~ ,fi) o('~ x ) and e^'® ' 



are added to "factor out" the free oscillation of the EM field, 
which has oscillation frequencies near (Oo. 
By imposing normalization conditions of 

m (0l(jC,Ol0lO«,O>ni = m(^!(Ol02(O)m=l, (10) 

the probability for finding the photon at location x (with x < 
indicating a photon propagating towards the cavity, and x > 
a photon propagating away from the cavity) is given by 

p Y (x,t) = \f(x,t)\ 2 (11) 

while the probability that the photon is in the cavity is given 
by |o£ 2 (f)|. In this way, the normalization condition of the 
joint quantum state 

|_J/(*,?)| 2 dx + |a 2 (f)| = l (12) 

is simply a statement about the conservation of total probabil- 
ity. 

The function f(x,t) can be viewed as the out-of-cavity pho- 
ton's wave function, while 1 0i (x, t)) m for each x can be viewed 
as the oscillator state that is entangled with each possibility 
for the out-of-cavity photon. On the other hand, a(t) can be 
viewed as the probability amplitude of the cavity mode, while 
|02)m can be viewed as the oscillator state that is entangled 
with the in-cavity photon. 

To facilitate calculation, for any joint quantum state 
we define 

I yi (*,*)>» = r (x\w)e iao{t - x) =f(x,t)\(l>i(x,t)) m (13) 
\yn(t)) m = {QH^e^' = a{t)\Ut)) m (14) 

Here \yi{x,t)) m , — 00 < x < is a series of vectors, 
parametrized by x, in the Hilbert space of the mechanical os- 
cillator, while I \j/2(x,t)) is a single vector in the Hilbert space 
of the mechanical oscillator. They together carry the full in- 
formation of the quantum state of the entire system. To further 
appreciate the role of |y/i) m and |y/2)m, we can project each 
of them into the position eigenstate of the oscillator, \y) m , ob- 
taining 

®l{t,x,y)= m (y\\l/i) m = f(x,t)<pi(y,x,t) (15) 
® 2 (t,y) = m(y\Y2)m = a(t)fa(y,t) (16) 

which can be viewed as the joint wave functions of the pro- 
jection of the entire state into J#i+ © M\- and J^f, respec- 
tively. Note that although f(x,t) and \$i(x,t)) m [and similarly 
a(t) and |02(f))m] share a phase ambiguity, \\j/i(x,t)) m and 
IV^(0)m> an( l nence < I ) i( f ! x i3 ; ) an d ^2(1, y) are well defined 
without ambiguity. 

C. Initial, Final States and Photodetection 

As special cases, we consider the quantum state of the sys- 
tem at t = (the initial state), and at very late times (the fi- 
nal state). For the initial state, the photon is propagating to- 
wards the cavity, and the cavity is empty. This corresponds 
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to a(0) = 0, and f(x,0) = 0. In particular, we also presume 
the initial state to be separable between the photon and the 
oscillator, with 

\W{0)) = J e iw ^F{x)\x) y dx®\^) m (17) 

Here F(x) is the slowly- varying part of the initial wave func- 
tion of the photon, and |0o)m the initial wave function of the 
oscillator. In other words, we have 

IVi(' = 0)) m = F(je)|to) m (18) 

|yi(f = o))m = o (19) 

with F(x) — for x > 0. At a sufficiently late time T, the 
photon will leave the cavity with unity probability, and we 
expect a(T) = and /(jc, T) = for x < 0. Mathematically, 

\Yi(x,t > T)) m = F out (x,t)\<j>(x,t)) m 
\Wi(t>T)) m = (20) 

with F out (t,x) = for x < and t > T. This is an explicitly 
entangled state between the out-going photon and the mirror, 
if | <j) (x, t ) ) m for different values of x are not all proportional to 
the same state vector. 

At an intermediate time t > 0, suppose a photodetector is 
placed at x = L > (i.e., for out-going photons from the cav- 
ity), then the probability density for photon arrival time at T 
is given by 

p L (T)= m {Yi{L,T)\Yi(L,T)) m . (21) 

In addition, by detecting a photon at this particular instant, the 
oscillator is left at a condition quantum state of \<j)(x,T)) m 

D. Evolution of the photon-mirror quantum state 

Applying the operations y(x\ and y(Q\a onto the (joint) 
Schrodingier equation 

ih d \ll = H\ V ) (22) 

we will obtain coupled equations for |y/i) m and |y/2) m . 
Throughout this section, we will mostly encounter states in 
the oscillator's Hilbert space, therefore we will ignore the sub- 
script "m" unless otherwise necessary. 

1. Free Evolution 

For by applying y(x\ to both sides of Eq. (22) we ob- 
tain 

[d t + d x + iH,„] \Wi(x,t)) = -y/fd(x)\\i/ 2 (t)). (23) 

Equation (23), without the 5-function term, simply describes 
the propagation of the initial photon towards the cavity, and 
the free evolution of the oscillator. This is because when the 
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Figure 3: (Color Online) Three regions of the t-x plane and the free 
evolutions of In region (i), the photon has not yet entered the 
cavity; the joint quantum state of the system is a simple free evolution 
of the initial quantum state, specified on / = 0, x < (green horizon- 
tal half line), see Eq. (25). In region (ii), the photon and the oscilla- 
tor evolve freely after propagates after the photon emerges from the 
cavity; the joint wave function depends on the wave function along 
x = 0, t > (green vertical half line). The red line diving regions (i) 
and (ii) corresponds to the 5-function in Eq. (23), which embodies 
the interaction between the outside photon and the in-cavity photon. 
Region (iii) is causally irrelevant to our experiment. 

single photon is outside the cavity, its propagation is free, 
while the oscillator's evolution is unaffected by light. 

Equation (23), is a first-order partial differential equation 
with characteristics along x — t = const. We hereby divide 
the t > region of the t-x plane into three regions: (i) x < 0, 
(ii) x > and t > x, and (iii) x > and t < x, as shown in 
Fig. 3. We can discard region (iii) right away, because it is 
not causally connected with our experiment. In the interiors 
of regions (i) and (ii) separately, Eq. (23) has the following 
general solution, 

\ Wl (x,t))=e-i»"^\C(t-x)) (24) 

with |C(v)) an arbitrary state-valued function of v. 

In region (i), |C(v)) can be specified by initial data along 
the half line of t — 0, x < 0; by using Eq. (24) twice, at (t,x) 
and (0,x — t), we obtain [See Fig. 3]: 

|V/i(x<0,f)) =F(jc-f)0 IB (f)|fo). (25) 

Here U m is the evolution operator for the free oscillator, given 
by 

G m (t)=e- i ^ t . (26) 
In terms of the Fock states \n), we have 

f/,„( f ) = £l"> e "' ( " + 'H>l (27) 

n 

Equation (25) corresponds to the photon's wave packet freely 
propagating along the positive direction of the x axis and the 
mechanical oscillator independently evolving under its own 
Hamiltonian. 

In region (ii), |C(v)) is specified by boundary data along the 
half line of x = 0+, t > 0, which we denote by 

\ Wl (t)) 0+ = \Y 1 (0+,t)). (28) 
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By using Eq. (24) twice, at (t ,x) and (t — x,0), we obtain 

\yr 1 {x>0,t)) = G m (x)\}fr 1 {t-x))^ (29) 

Henceforth in the paper, 0+ and 0— stand for x H> 0+ (x ap- 
proaches from positive side of the axis) and x n- 0— (x ap- 
proaches from negative side of the axis) respectively. Equa- 
tion (29) corresponds to the free evolution of the out-going 
photon and the mechanical oscillator. 

2. Junction Condition 

The 5-function on the right-hand side of Eq. (23) relates the 
out-going photon to the decay of the in-cavity photon and the 
reflection of the in-going photon. To take this into account, we 
simply integrate both sides from x = 0— to x = 0+, obtaining: 

|^i(0+,0) = |^i(0-0)-V7l^2(0) (30) 

This expresses the out-going wave as a combination of the 
promptly reflected incoming wave and the wave coming out 
from the cavity. 

3. Coupled Evolution 



which indicates an overall down-shift of eigenfrequency. Here 
we have further defined the displacement operator 

fi(j3) = exp[j3(i t -b)} (36) 

with b and b^ the annihilation and creation operators for the 
free mechanical oscillator (i.e., before it couples to light). As 
we shall see in Sec. VI, /3 will become an important charac- 
terizing parameter of our optomechanical device; for example, 
/3 > 1 is the regime in which the device is nonlinear. 

For the photon, Eq. (31) means that the in-cavity photon 
is continuously driven by the in-coming photon (right-hand 
side) and decays towards the out-going photon (as indicated 
by the y/2 term in the bracket on the left-hand side). The 
above discussion, together with the initial data of \\jf2) = at 
t = gives 

\Vi) = VY fe-i^U r (t-t>)\ Wl (t>)) _ (37) 
Jo 

where 

Uy(t) = = £ \h)e-^ n+l l 2 -^ m "''(hl (38) 

n 

which is the modified evolution operator of the oscillator 
when the photon is in the cavity. 



By applying Y (Q\a to both sides of Eq. (22) and using 
Eq. (30), we obtain: 



d, + ^ + iH r 



(31) 



Here as in Eq. (30), we have defined |yi)o± = |Vi(0±,f)). 
We have also defined 



A _ Py , mcoj(y-a) 2 2 
Hr= 2^ + 2 P W " 



with 



a 



m(£>; 



CO m ^/2m(0,„ 



(32) 



(33) 



The operator Hy can be viewed as the modified Hamiltonian 
for the mirror when the photon is present in the cavity. Here a 
characterizes the shift in equilibrium position of the harmonic 
oscillator when the photon is inside the cavity and applies a 
constant force to the oscillator, while /3 (as seen from this 
equation) modifies the eigenfrequency of the harmonic oscil- 
lator. It is easy to work out the eigenstates and eigenvalues of 
Hf. the eigenstates are 



\n) =e ia Py\n) =£(j3)|n) 



(34) 



which are simply displaced from the original Fock states in 
phase space, due to the change of equilibrium position, with 



H Y \n)= {n+ 1 --j3 : 



(O m \n) 



(35) 



4. Full Evolution 

The full evolution of the entire system's quantum state can 
now be obtained by combining Eqs. (29), (25), (30) and (37). 
In order to study the out-going photon, we only need to con- 
sider the region x > and t > x (see Fig. 3), because the 
emerges from the cavity at t > 0, and it propagates with c = 1 . 
For this region, we obtain a compact-form solution of 

\Wi(x,t)) = M\<h) (39) 
where |0o) is the initial quantum state of the oscillator, and 



Jl/= J o g{t-x,t')U m (x)U Y (t-x-t')tJ m (t')dt' 

J\ iP 2 co m (t-x-t') 







where 



with 



g(t-x,t')e l i 
D(pe- ,a "' x )D{-lie' a "' {, '-^)U m (t)c\t' . (40) 

g (t,t') = G(t-t')F(-t') (41) 



G(t) = 8 + (t) + ye- 



(42) 



the cavity's optical Green function. Here the subscript + for 
the 5 function indicates that it's support lies completely in 
the region t > 0. Within the operator M [Eq. (40)], the factor 
g contains two terms, the first contains a 5-function and the 
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second an exponential decay over time. The first term corre- 
sponds to the photon being promptly reflected by the cavity's 
front mirror, while the second term corresponds to the pho- 
ton staying inside the cavity, for an amount of time equal to 
t — x — t', which ranges from to t — x. As a sanity check, 
it is straightforward to see that when mass of the oscillator 
approaches infinity, Uy coincides with U m , and M simply de- 
scribes the photon's propagation and the independent evolu- 
tion of the oscillator. 



III. SINGLE-PHOTON INTERFEROMETER: VISIBILITY 

In this section, we will use the results of the previous sec- 
tion to analyze the single-photon interferometer. 

A. The configuration 

We consider a scheme proposed and analyzed by Marshall 
et al. [7], which is shown in Fig. 1. This Michelson interfer- 
ometer (with 50/50 beamsplitter) has two arms: in the north 
arm, the end mirror in cavity A is movable, and initially pre- 
pared at a quantum state |</>o), whereas mirrors in cavity B or 
east arm are fixed. We assume the photon is injected from the 
west port, while a fixed photodetector is placed at the south 
port. Apart from mirror A being movable, the two cavities are 
otherwise identical: with the same input-mirror power trans- 
missivity T, length L (for cavity A, counted from the zero- 
point of A's displacement). The front mirrors are placed at 
equal macroscopic distance from the beamsplitter, while there 
is a phase detuning of <p in arm B for C0q 1 . In our convention, 
if mirror A is at zero point and <p = 0, the photon will always 
return to the west port. Henceforth in the paper, we shall refer 
to the west port as the input port, and the south port the out- 
put port — although we may not always find the photon at the 
output port. Indeed, whether and when the photon arrives at 
the photodetector is jointly determined by <p and the state of 
motion of mirror A. 

In particular, we shall use p(t) to denote the probability 
density for the photon to arrive at the detector at t (which can 
be measured by repeating the experiment many times). If we 
idealize the arrival time of the in-going photon (at the front 
mirror) to be t = 0, and ignore the macroscopic distance be- 
tween the front mirrors, the beamsplitter, and the photodetec- 
tor, then we are interested in p(t) at t > 0. We further define 
an instantaneous fringe visibility 

V{j) = p^(t)-P^(?) 

Pmax(f) +Pmin(0 

which measures the degree of coherence between the two 
components of returning photons at the beamsplitter, and can 



To give rise to a detuning, we assume that all optical frequencies we con- 
sider are centered around (Bo, and we offset the location of cavity B from 
symmetry by a length / such that (OqI = (p/2. 



only become unity if at time t the joint mirror-photon quantum 
state is separable, as we shall see more clearly in Sec. HID. 

B. The role of the beamsplitter and a decomposition of field 
degrees of freedom 

In Sec. IIC, we have studied in detail how the photon first 
affects the x < components of the optical field out-side of a 
cavity, then interacts with the mirror, and finally returns back 
to the x > components of the optical field. The scenario 
for a Michelson interferometer is slightly more complicated: 
we now need to consider a set of input fields that replaces the 
x < single field in the single-cavity case, and a set of output 
fields which replaces the x > single field. 

As shown in Fig. 4, the annihilation operators of the input 
field for the two cavities are (j_, a_, b-), while those of 
the output fields for the cavities are k + , d + , b+). Each 
of these files are defined as a function of — °° < x < +°°, with 
x = corresponding to the position of the beamsplitter, and 
positive direction along the arrow shown in Fig. 4. Ultimately, 
we need to calculate the fields of j + and k + in terms of y_ and 

Note that at by allowing x to run through the entire real axis, 
we have assigned two input fields and two output fields to each 
point along the optical path (note here that "input" and "out- 
put" refer to the cavities, not the beamsplitter). This redun- 
dancy is necessary for a simplified treatment of the beamsplit- 
ter: instead of treating its internal dynamics, we simply view 
it as a mapping between the two different representations of 
the input and output fields. One representation (j±, k±) corre- 
sponds to the point of view of observers at the west and south 
ports, pretending that the beamsplitter does not exist; the other 
(a±, b±) corresponds to the point of view of observers at the 
east and north ports. 

The conversion between the two representations takes the 
same form as the "input-output relation" of the beamsplitter: 

? / \ ]±{x)-k ± (x) } ± {x)+k±(x) 

= — — ' fl ± w = — — • ( 44 ) 

As an example, consider a quantum state in which a (instanta- 
neous) photon is injected from the input port, which, accord- 
ing to the mapping in Eq. (44), has two equivalent representa- 
tions: 

<■+ , m , a* (xn) (xc\) , , 

£(*o)|0) = V i°) • (45) 

As time grows, the quantum state evolves as xq — > xq + t. At 
any instant, the left-hand side represents a single photon prop- 
agating from west to east, and continue through the location of 
the beamsplitter. The right-hand side represents a photon has 
a two-component wave function, the first component propa- 
gates northwards, the second eastwards. 

Although the two representations are equivalent, we still 
prefer to use the south-west representation when treating the 
generation and detection of photons, and the north-east rep- 
resentation when treating the light's interaction with the cavi- 
ties. 
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J+ 



Figure 4: We illustrate the fields entering and exiting each of the four 
ports of the interferometer. We use arrows to define the positive sense 
of the coordinate used to label their locations. For each of them x = 
corresponds to the location of the beamsplitter. 



C. Interactions between light and cavities 

For each individual cavity, we intend to apply the result of 
Sec. II A. We note that a_ (x) (for x < 0) and a + (x) (for x > 0) 
defined in this section maps to the c(x) (for x < 0) and d(x) 
(for x > 0) respectively, as defined in Sec. II A and illustrated 
in Fig. 2. For this reason, we define 



and 



a(x) = 



b(x) 



d-(x), x < , 




(46) 



(47) 



b+[x), x> 0. 



In this way, a(x) and b(x) here both map to c(x) defined in 
Sec. II A. [The a and b here are not to be confused with op- 
erators of the optical mode and the mechanical oscillator — 
we shall always explicitly include the argument (x) for these 
continuum operators.] We further define 



and 



/(*) = 



k(x) = 



J-(x), x<0, 

/+(*), x>0, 

k-(x), x < , 

k+ (x) , x > . 



(48) 



(49) 



Furthermore, for fields a, b, j and k, the transformation rela- 
tions Ea. (44) also apply. 

Now suppose at t = 0, we have a photon coming from the 
input (west) port with arbitrary wave function F(x) [like in 
Eq. (17), here F(x) =0 for x > 0]. The initial quantum state 
of the entire optomechanical system is 



|V/(0)) = / d*F(jt)/L(*)|0) y ®|to>, 



(50) 



Since we would like to investigate this state's evolution when 
the photon reaches the cavity, we covert into the north-east 
representation: 



Mo)) = ^[|v'A(o)> + |vfe(o)>] 



(51) 



Here we have defined 



IVa(O)) 



lv*(o)> 



dxf(x)fl t (^)|O) r ®|0 o > A (52) 
dxF(x)b f {x)\Q) Y ®\<l>o) A (53) 



in which we have already taken Eqs. (46) and (47) into ac- 
count. 

Here |i//^(0)) corresponds to the case in which the pho- 
ton enters Cavity A with the movable mirror, and | i//s(0)) the 
case in which the photon enters Cavity B with the fixed mir- 
ror. As time goes on, these two states evolve individually, and 
Eq. (51) remains true for t > 0. For the Cavity A component 
| V)a, we have [cf. Eq. (39)] 

|Vk(0)= I' X dt'g(t -x,t')e^ 2co '" {t - x -^D(pe- ico '" x ) 
Jo 

D{-^e iam( f'-^)&\x)\^{t)) (54) 
where we have defined 

\^{t)) = U m {t)\0) y \^) A , (55) 
while for | y/) b, we set /3 — > and obtain 

|Vfe(r)) =e i * r X 6t'g(t-x,t')P(x)Mt)). (56) 
Jo 

D. The final state 

In order to describe the quantum state seen by the photode- 
tector, we map a and b into j and k, only keeping the k com- 
ponent. We further project onto the single-photon basis of 
y(0\k(x), assuming x = 0+, obtaining 

\V(t))m=\ [|V&(0>m + ^V&(0>m] (57) 



with 



\WA(t)) m = / Al'g{t,t')0{t-t')\<fo{t)) (58) 



with 



and 



Wn{t)) m = / dt'g(t,t')\<j> (t)) 

Jo 



\fa{t)) = O m (t)\fa) 



6{t) = e^^'Di^Di-Pe-" "-') , 



(59) 



(60) 



(61) 
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in particular 0(0) = 1. In this way, we are using the same 
notation as Eq. (39), and we can use Sec. II C for obtaining 
photo-detection probability density at each time t > 0, which 
is given by 



Pit) 



Wa) 



|i/A B ) m || 2 + 2Re (^(VaIVb)™) 



(62) 



which, when adjusting values of <p, leads to an instantaneous 
visibility of [Cf. (43)]: 



v(t) = 



2|m(VA|Vfi).: 



\\¥A)r 



\\Vb)v 



12 ' 



(63) 



It relies on how different is from ys, which indicates how 
much the movable mirror in Cavity A is capable of "learning" 
about the existence of the photon in Cavity A. At any instant, 
if y/4 is proportional to y/g (differ by a phase), the state of the 
movable mirror does not change, and therefore we have a per- 
fect visibility. By contrast, if the photon is able to transform 
the movable mirror into a state substantially different from its 
freely evolving state, e.g., the orthogonal state in the extreme 
case, then we will have a significantly reduced visibility. 

Similar to Eq. (39), here y//t and y/s each has a promptly re- 
flected part [which arises from the 5-function part of g(t,t')], 
and a part in which the photon enters the cavity [which arises 
from the exponential decay part of g(t,t')]. It is the second 
part that contributes to the reduction of visibility. 



E. Examples 

We consider an experimental situation with the central fre- 
quency of the injecting photon tuned to the resonant frequency 
of the cavity, with a wave function of 



F(x) = V2Te Tx ®{-x). 



(64) 



Here F measures the frequency-domain width of the photon. 
We further assume that the mechanical oscillator's eigenfre- 
quency (when uncoupled with light) is equal to the cavity 
bandwidth, or co m = y. As in Ref. [7], we assume that the 
mechanical oscillator, i.e., the mirror, is initially prepared at 
its ground state: 



\<h) = \o)a. 

With these specializations, we have 

\YA(t)) m =C(t)[\0) + Y\M(t))} 



|VCi(0)m = C (0 



|0) + y/V/(*-0|0) 
Jo 



with 



C(f) = V2Te- (r+i 



m m /2)t 



,<X-y/2)t 



(65) 

(66) 
(67) 

(68) 



and 



\M(t)) = f fit-tyiPWit-^-smcOnit-t 1 )] 

Jo 



(69) 



By comparing with Sec. HID, we first find that visibility de- 
pends on the similarity between \M(t)) and its counterpart in 
Eq. (67): when they are similar to each other (e.g., when j3 < 1 
) or when they do not contribute significantly to \y>A.B{t))m, 
the visibility will tend to be high. By contrast, in order to 
achieve a complete incoherence, we need \M(t)) to contribute 
significantly, and nearly orthogonal to |0) — and this requires 
/3 > 1 . The arrival probability density (62) and contrast defect 
(63) can be computed if we use 



<0|j3> = <0|D(/3)|0>=* 



(70) 



In Fig. 5, we plot maximum and minimum of the probabil- 
ity density in the left panels, and visibility in the right panels, 
both as functions of time. We have chosen j3 = 0.5 for upper 
panels, /3 = 1 .2 for middle panels and j3 = 2 for lower panels. 
In each panel, we have also shown curves with F = 0.2 (red 
dotted), r = 1 (blue dashed) and T = 2 (solid black). As j3 in- 
creases (as we move from upper to lower panels), the photon's 
ponderomotive effect on the movable mirror increases, there- 
fore the visibility is able to vary more. This means /3 > 1 is 
necessary (but not sufficient, see below) for visibility to sub- 
stantially decay and then revive — a feature Ref. [7] has used 
to search for decoherence effects. 

On the other hand, another condition for visibility to first 
decrease and then revive, and repeat on, seems to be F > 1, 
as also indicated by each of the right panels of Fig. 5. In 
addition, as F ^> 1 , our result becomes comparable to Ref. [7] . 
Qualitatively, this is because for F 3> 1, if photon does arrive 
at a time around ~ 1, we can be sure the photon has interacted 
with the mirror — and we can roughly treat the photon as 
already within the cavity at t = 0. 

Mathematically, for t 3> 1 /F, the conditional quantum state 
of the mirror given photon detection at time t could be approx- 
imately written as : 



IV)m = 



|_ e -(y+i0i*)f/2 

2a 

■ i( P\0) +e^ 2mm 'D(j5)\-j5e- i0> »' t ) 
y_ e -(y+ia m )t/2 



2a 



?" P |0) -^g'f'^nt-s.mamt) \p _p e -i(0„,t\ 



(71) 



(72) 



This is consistent with results of Ref. [7]. 

However, in order for a ^ 1 and to observe a revival of 
visibility, we have to wait till t > 2n. The probability for de- 
tecting the photon at such late times is exponentially small — 
as indicated by the left panels of Fig. 5 This means we may 
have to make a trade off between having a very sharp revival 
of visibility and being robust against loss and able to cumulate 
enough statistics within reasonable amount of time. 



IV. CONDITIONAL QUANTUM-STATE PREPARATION 

In this section, we show how to engineer an arbitrary quan- 
tum state of the mechanical oscillator by injecting a single 
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Figure 5: (color online) (left) Probability density and (right) fringe 
visibility for the photon to come out with different /3: (Top-to- 
bottom: first row, j3 = 0.5; second row, /3 = 1.2; third row, jS = 2). 
For each jS, three different values of T are considered for compar- 
ison: r = 0.2 (red dotted), l(blue dashed), 2(black solid). All the 
calculation assume 7=1, (O m = 1 . For probability density plot, the 
upper line of the same color is the maximum value of the probability 
density, the lower one is the minimum value. 



photon with specifically designed wave function and by post 
selecting the arrival time of the output photon. Note that un- 
like Refs. [10, 11], our state preparation procedure is condi- 
tional. This guarantees a pure quantum state for the mechan- 
ical oscillator, but requires a low decoherence rate and a high 
detection quantum efficiency for the out-going photon. 



A 



Single-Photon 
Source 



B 



Detection 



Figure 6: (color online) The sample device which uses single photon 
to prepare mechanical oscillator quantum state. Here the the detuning 
phase for the mirror on the east arm is adjusted such that the promptly 
reflected photon will come out from west port, with probability 
coming out from south port. 



replaced by 8{t — t'), as we have a perfectly reflecting mirror 
instead of a cavity here, namely, 



|Vu(O)m = |0D(O> 



(73) 



To proceed, we further adjust the detuning phase <p in 
Eq. (57) such that at the dark port, the promptly reflected wave 
from the front mirror of cavity A exactly cancels the promptly 
reflected wave from the mirror B. In this case, having a pho- 
ton emerging from our detection port (Fig. 6) automatically 
indicates that the photon has entered the cavity and interacted 
with the mirror; Eq. (57) or the conditional quantum state of 
the mechanical oscillator (unnormalized) is given by: 



1 



with 



Mt)) m = - / dt' gp (t,t')o(t-t')\fo(t')) 



gp (t,t')= Y e-^ t -^F(-t'). 



(74) 



(75) 



As g p (t,t') is related to the input photon wave function F(x), 
by modifying input photon wave function, we can therefore 
engineer the conditioning mechanical oscillator quantum state 
\y(t)) m . Even if there is a finite probability that the photon 
will come out through the west arm or the bright port, once 
we detect a photon at time t at the dark port, we know that it 
must come from arm A and it also has stayed in the cavity A 
for a certain amount of time. 



A. The configuration 

The scheme is shown in Fig. 6. It is very similar to the 
single-photon interferometer discussed in the previous sec- 
tion, except that in the east arm we replace the cavity B with 
a perfectly reflected mirror B. In this case, most of the pre- 
vious analysis are still valid: Eq. (57) to Eq. (59). The only 
difference is that the g(t,t') function in Eq. (59) needs to be 



B. Preparation of a single displaced-Fock state 

First of all, we notice that when different in-coming photon 
wave function F's are used, if we keep conditioning over the 
same photon arrival time t , the conditional quantum states we 
obtain for the mechanical oscillator will depend linearly on 
F. In other words, if F\ allows us to prepare and F^ 
allows us to prepare \<j>2), then injecting a new photon with a 
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Figure 7: A sketch of the phase-space trajectory of the mechani- 
cal oscillator. The Wigner function of the initial state |0) is repre- 
sented by the shaded disk, the dot marked with on the real axis is 
the new equilibrium position of the oscillator when the photon is in 
the cavity, while the dashed circle is the trajectory of the oscillator's 
Wigner function when the photon is inside the cavity. Detection of 
the out-going photon at t = 2nlt/a> m corresponds to superimposing 
all mechanical-oscillator quantum states along the dashed trajectory, 
weighted by the photon's wave function. 



\ V (t)) m = \n)=£>(p)\n) 



0,1,2. 



(76) 



These displaced Fock states are simply Fock states of the os- 
cillator when the photon is inside the cavity, see Eq. (34). 

Let us assume that the mechanical oscillator is initially pre- 
pared at its ground state. Before studying preparation of an 
arbitrary conditional quantum state for the mechanical oscil- 
lator, we first show that we can prepare a conditional state with 
an arbitrary quantum number n, by injecting a photon with the 
following wave function: 



(r/2-ij3 2 a> m +ma>„, 



]x &{-x). 



(77) 



As we plug Eq.(77) into Eq.(74) we obtain the conditional 
quantum state of 



\W{t))r, 



D(P)f/ 2 e-i t+i P 

2(O m 



dT'e- inx '\-pe' { 



(78) 

with T = cot. This is a coherent superposition coherent states, 
which in the complex amplitude domain all line up in a circle 
with radius j3 around the center located at complex amplitude 
equal to j3; these states are parametrized mathematically by 
D(f})\ — fie"!'). These states are superposed with the same 
magnitude, but different phases, due to the decay rate of y/2 
in the F chosen by Eq. (77). Obtaining such a state is under- 
standable, as given the photon detection at t, the actually time 



t' for the photon staying inside the cavity is uncertain, and we 
have to sum up all the possible contributions from to t . This 
situation is illustrated in Fig. 7. 

One important feature in the above expression is that the 
integrand is a periodic function. If we denote 



t = co,„t = 2nN + A(j) , 



(79) 



where N is some integer and A(j) is the residual phase ranging 
from to 2tz. In this way, the integral in Eq. (78) then becomes 



rln 




N / 4 




Jo 


Jo 



d(j)t 



(80) 



In the limit of N 3> 1, when the photon arrives at the pho- 
todetector with a delay large compared to the oscillator's os- 
cillation period, the first term in Eq. (80) always dominates. 
This means we obtain the same conditional state if we restrict 
T around an integer multiple of 2n, or make sure it is large 
enough. This leads to the interesting effect that in the asymp- 
totic limit of T — > +°°, the conditional state will be indepen- 
dent from T. In practice, however, although the integral (80) 
increases with N, the exponential decay factor in Eq. (78) al- 
ways favors simply choosing N = 1. It is straightforward to 
evaluate this conditional state; using 



i-2n 
Jo ' 



— in& „e v a ] 



l 



superimposed wavefunction F = (XiFi + CX2F2 will allow us to we have 
prepare «i |<^i) + «2|<fc} 

This means we only need to show how members of a com- 
plete basis can be prepared, and we choose this to be 



Ik 



d</>< 



■pe* 



\0)=-W)*\0), 



2n{-Pfe- P -r 
2K\n){n\-p), 



(81) 



which means 



IV)r 



<>>m ( 



(82) 



(83) 



This is indeed proportional to \h), as promised. Here we have 
used 



-PI") 



(84) 



Since the probability for the returning photon to arrive at 
precisely 2%jco m is zero, we must allow an interval around 
this target, which on the one hand provides us with a non-zero 
probability, but on the other hand makes the conditional state 
imprecise. If we require the actual conditional state to have a 
high overlap with the target state (or high fidelity), 



> 1 



(85) 



then, by perturbing the integration upper bound of Eq. (82), 
we obtain the following requirement on the allowed photon 
arrival time 



h-2n\ < At = V8n 2 e 



-P\n) 



(86) 
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Figure 8: Probability for obtaining displaced Fock states |1) (red 
solid), 1 2) (blue dashed), |5) (magenta dotted) and 1 10) (black dash- 
dotted), a range of j6 and minimum state overlap of 1 — e. Vertical 
gridlines are draw for /3 = 1 , \/2, v5 and \f\0; these are the locations 
where maxima of -Pi^.s.io are reached. 



which, for each trial of the experiment, would happen with a 
probability of 



P = |m(vlV> 

= 2V8e 



|2 2At 



\(-P\") 



(87) 



And this would be the probability with which we can create a 
conditional state with overlap at least 1 — e with the target. 
From Eq. (87), we further notice that we should fix 



y/to m = 3/{2n) 
in order to obtain a maximized success probability of 

27 |(-j3|«)| 3 



P n = V8e 



4*Vl-|H3|»> 



(88) 



(89) 



For each «, the maximum of P„ is reached at j3 = y/h. In Fig. 9, 
we plot P n for a range of j5, for e = 0.1, or a state overlap 
of > 90%. We can see that the probability of producing \h) 
decreases rather quickly as n increases. 

This dependence (89) on j3 comes from two sources, which 
we can understand better by going to the phase-space refer- 
ence frame centered at the equilibrium position of the oscil- 
lator when the photon is inside the cavity. In this reference 
frame, the complex amplitude of the coherent states being su- 
perimposed are located on a circle with distance j3 away from 
the center, while the target we would like to prepare is sim- 
ply the Fock state \n). Although the photon's wave function 
selects out an oscillator state proportional to \n), this post- 
selection does not improve the intrinsic overlap between all 



those that participate the superposition, which is actually pro- 
portional to 



\(-P* 



= \(-P\n) 



(90) 



This explains the dependence of m (y| y/) m on beta. The other 
factor of dependence on j3 is that when the target state has a 
very low overlap with the individual members \^e'^) of the 
superposition, the requirement on the accuracy of photon ar- 
rival time, or At, increases, as shown in Eq. (86). 



C. Preparation of an arbitrary state 

Since the displaced number states forms a complete basis 
we can expand any target state as 



n=0 



n=0 



(91) 



Since a linear combination of F's leads to a linear combination 
of conditional states, we simply need to apply the result of the 
last subsection and have 



/y e {r/2-iP 2 (O m )x +oo 

F ( X ) = ~ y E ^ 

^ n=0 



(92) 



with 



1 1/2 



\-P\n) 



(93) 
(94) 



This is an additional periodic modulation (with period 
27r/o) m ) of the photon's wave function. We caution that in 
order for the summation in Eq. (92) to converge, if c„ does not 
go to zero for all n > N, then it must decay very fast when 
n — > +°o, due to the presence of the \fn\ factor (which grows 
faster than j3~"). 

As in the previous subsection, we obtain the conditional 
state at T = CO m t — 2%, Alt, as well as any T that is sub- 
stantially large. Again, let us consider T = 2n, this gives the 
conditional state of 



IV)m = 



co,„Z 



IVtg) 



(95) 



We can use the same approach as the previous subsection 
to evaluate the probability with which this conditional state is 
achieved with a high overlap. For a minimum overlap of 1 — e, 
we require 



\2n-x\ < At = 



L^r/mlv'l-IHSIV'tg}! 2 



(96) 



Note that this At diverges if L,„~ c m = 0, because in this case 
the overlap does not vary at 0\(x— 2n) 2 ] order. Assuming the 
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Figure 9: Minimum success probability for states in Hilbert spaces Mip. ...7 (solid curves with markers), together with success probability 
for producing single displaced Fock states, 1 2 ... 7 (dashed curves without markers). Fidelity is fixed at 10%. Note that Pq would become 
greater than 1 at low values of j3 — but in this case our approximation in obtaining At breaks down. 



target state to be generic, then the probability for obtaining 
this state is then 



P\ w) =2V8e 



lity - 



2* I 2 



i-|i„ro(-j3i") 2 ^ 



-1/2 



J k 



-""^tol + i^ 



(97) 

Here the choice of 7/ co m depends on the target quantum state, 
but if we assume this dependence is weaker than the pre- 
factor, and continue to use Eq. (88), then we obtain 



27 

73 



i-nro(-^i") 2 ^ 



2 = I 2 



-1/2 



(98) 



-m=0 L '" I 



' j-k=0 1 + 



iKi(j-k) 



As it turns out, P| ^ depends on the detail of 1 1//) — even if 
we only try to create a combination of |0) and |1), the combi- 
nation coefficients would lead to very different success prob- 
abilities. In order to provide a concrete measure of the ability 
of our state-preparation scheme, we have chosen to compute 
the minimum success probabilities of creating all the states 
in the mechanical oscillator's Hilbert subspaces spanned by 
the lowest displaced Fock states, e.g., J#\ = Sp{|0),|l)}, 
J% = Sp{|0), |1), |2)}, etc. We define 



Fock states. For smaller /3, it is much more difficult to pre- 
pare a higher Fock state, therefore, if /3 is sufficiently small, 
the difficulty of preparing Jtfj is dominated by the preparation 
of I J), the single most difficult state in the space to prepare — 
and therefore Pjf. agrees with Pj. Vice versa, for sufficiently 
large j3, the difficulty of preparing J(?j lies in the preparation 
of 0), and therefore P^. would agree with Pq. The fast tran- 
sition between the two extremes indicates that when trying to 
prepare states in Jtfj, the difficulty either lies in |0), or in | J), 
and only for a small region of j3 the two difficulties might 
compete with each other — while none of the intermediate 
states contribute to the difficulty of state preparation. This is 
consistent with the relative locations of the P n curves in Fig. 9: 
(i) for any /3, P\ 2,...j-i are always much greater than the min- 
imum of Pq and Pj, and (ii) as we move away from the /3 
at which Po and Pj crosses each other, their discrepancy in- 
creases quickly. 

As a matter of practicality, we see that if we choose j3 ss 
0.87 the probability of achieving, with an overlap (or fidelity) 
above 90%, any superposition of |0) and 1 1) (i.e., any member 
of the subspace is guaranteed to be above 6.3%. On the 
other end, with a probability of at least 0.1%, we can produce 
all states in the 8-dimensional subspace of ^7. 



V. PRACTICAL CONSIDERATIONS 



min Pi 

\y)ejfj 



¥i ■ 



,/=o 



a, G 1 



(99) 



In Fig. 9, we plot Pj&, P^ 2 , P.^ 7 as functions of (in 
solid purple curves). Because C ^2 C • ■ • C ^7, it is in- 
creasingly difficult to create all states in Jtffj with higher values 
of /', and therefore P^ > P.^ 2 > . . .Pj&, namely our success 
probability decreases globally when j increases. In fact, as we 
overlay the single-Fock-state success probabilities Pq,P\,..., 
P5, we also discover that for any P^(j3), it asymptotes to Po 
at higher j3, and to Pj at lower j3; moreover, the transition be- 
tween these two asymptotic regions are brief, and the P^. ()3 ) 
curves do not lie much below the minimum of Po and Pj. 

This asymptotic behavior can be understood from the be- 
havior of P„, the success probability for single (displaced) 



In order to realize such a state-preparation scheme, we need 
to fulfill the following three requirements. The first require- 
ment is that the series in Eq. (92) be converging. This can be 
satisfied if J3 > 1. To see what this means, we restore all the 
physical units: 



k/(2(0,„) 
^hmco m /2 



hcoQ 



hm(O m 



2o)„, 



(100) 



It characterizes the momentum kick of photon hcoo/c to 
the oscillator during one oscillation period compared to the 
ground state momentum uncertainty ^hmO) m /2. The momen- 
tum kick from the photon needs to be big enough to substan- 
tially change the mirror state. The second requirement is that 
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cavity bandwidth be smaller than the mechanical frequency 

Y <C0, n . (101) 

This is because we need to wait at least several oscillation pe- 
riods to approach the asymptotic state, and the photon should 
be long enough such that we have a finite probability for de- 
tecting photon at t > <a~ 1 . 

Combining the above two conditions, we obtain the follow- 
ing relation 



9 



< 



n 



2m(O m 



(102) 



where A is the optical wavelength of the photon, & is the 
cavity finesse. This means the cavity linear dynamical range 
much be less than the zero point uncertainty to realize the op- 
tomechanical nonlinearity. An alternative scheme has been 
proposed to make it more achievable experimentally [16]. 

The third requirement is that the thermal decoherence effect 
be small within one mechanical oscillation period, namely [cf. 
also Eq. (5) in Ref. [7]]: 



Q> 



kT E 
hco,„ 



(103) 



where Q is the mechanical quality factor of the oscillator and 
Te is the environmental temperature. These three require- 
ments can be achieved experimentally, e.g., the current setups 
shown in Refs. [3, 4] the one proposed in Ref. [6]. 

Finally, we require the the capability of generating a single 
photon with an arbitrary wave function with duration compa- 
rable to the mechanical oscillation frequency of the photon. 
This is possible with cavity QED systems, as has been dis- 
cussed by Ref. [17-19]. 



decaying profile and with the movable mirror initially pre- 
pared at the ground state. This scheme has been proposed by 
Ref. [7] to explore decoherence of a macroscopic oscillator, 
although in that proposal the photon has been assumed to start 
off from inside the cavity. In the limit when the photon pulse is 
short (or a 3> y), we did recover the result of Ref. [7], although 
our result deviates significantly when a becomes comparable 
to y. We believe this is experimentally relevant, because in 
the case a> 7, the probability for the photon to exit from the 
detection port is very small, and therefore the experiment may 
suffer significantly from imperfections. 

We have also studied the use of such nonlinear optomechan- 
ical interactions to prepare the mechanical oscillator into an 
arbitrary quantum state — similar to the proposal of Ref. [12], 
although not having to require that the photon to start off from 
within the cavity. To realize this, we require that: (i) the 
optomechanical cavity must be working in the nonlinearity 
regime [i.e., the cavity's spatial line width must be less than 
the oscillator's zero-point position fluctuation, see discussions 
above Eq. (102)], (ii) the cavity's frequency width must be 
less than the mechanical oscillator's angular frequency, (iii) 
the thermal decoherence time must be less than several times 
the mirror's period of oscillation, and (iv) we must be able to 
engineer the single-photon wave function arbitrarily, at a time 
scale comparable to the mirror's oscillation period and with 
coherence time longer than the cavity storage time. Although 
we have shown mathematically that all quantum states whose 
expansion coefficients in the displaced Fock states \n) drop 
sufficiently fast as n — » +°° can be prepared by modulating 
the wave function of the incoming photon and conditioning 
over the arrival time of the returning photon, in practice we 
will be confined to the superposition of a handful of nearby 
displaced Fock states. 



VI. CONCLUSIONS 

We have presented an exact solution to the open quantum 
dynamics of an single-photon interferometer with a movable 
mirror. Since the photon number is preserved, we have been 
able to write the total wave function of photon as three compo- 
nents: incoming photon, inside-cavity photon and out-going 
photon. We analyzed the details of how the photon exchanges 
between the cavity mode and the external continuous field. 

We studied the fringe visibility of the interferometer in a 
specific case by injecting a single photon with exponentially 
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